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HIGHLIGHTS 


•  The  reaction  orders  of  V02+  and  H+  for  the  oxidation  reaction  of  V02+  are  calculated. 

•  A  new  reaction  mechanism  is  proposed  to  describe  the  oxidation  of  V02+. 

•  A  real  kinetic  equation  is  built  up  for  the  oxidation  reaction  of  V02+. 
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With  the  consideration  of  optimizing  the  performance  of  the  all-vanadium  redox  flow  battery  (VRB),  the 
oxidation  reaction  mechanism  of  V02+  on  a  rotating  graphite  disk  electrode  has  been  investigated  by 
potentiodynamic  polarization  in  sulfuric  acid  solutions  with  various  pH  and  vanadium  concentrations. 
Furthermore,  the  reaction  orders  of  V02+  and  H+  for  the  oxidation  reaction  of  V02+  have  been  calculated 
from  the  polarization  results  and  compared  with  the  theoretical  results  according  to  the  possible  reaction 
mechanisms  available  in  the  literature.  However,  a  new  oxidation  reaction  mechanism  has  been  pro¬ 
posed  to  describe  the  oxidation  of  V02+  at  last,  and  the  theoretic  reaction  orders  of  V02+  and  H+  based 
on  the  new  mechanism  are  consistent  with  the  experimental  results  when  the  electrochemical  reaction 
is  the  rate-limited  process.  Moreover,  a  corresponding  kinetic  equation  has  been  established  for  the 
oxidation  reaction  of  V02+  on  a  spectroscopically  pure  graphite  electrode,  and  can  be  well  used  to  predict 
the  polarization  behavior  in  V  (IV)  acidic  solutions. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

With  the  rapid  and  continuous  exhaustion  of  the  available  fossil 
fuel,  the  renewable  energy  resources  such  as  wind,  solar  and  tidal 
power  are  required  to  be  efficiently  developed  and  utilized  in 
future.  Due  to  the  random  and  intermittent  nature  of  these 
renewable  energy  resources,  energy  storage  technologies  should  be 
obtained  for  their  stabilization  and  smooth  output  [1].  To  date,  as 
an  electrochemical  energy  storage  system  the  all-vanadium  redox 
flow  battery  (VRB)  using  the  same  element  in  both  half-cells  has 
shown  a  great  promise  in  various  pre-  commercial  to  commercial 
stationary  applications  [2],  with  some  attractive  features  like  quick 
charging-discharging,  simple  structure,  flexible  design,  high  effi¬ 
ciency,  high  safety  and  long  cycle  life  [3].  In  the  VRB,  the  V0f/V02+ 
and  V2+/V3+  couples  in  concentrated  sulfuric  acid  solutions  are 
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applied  as  positive  and  negative  electrolytes,  respectively.  VRB  re¬ 
leases  or  stores  the  electric  energy  through  the  valence-state 
changes  of  vanadium  ions,  involving  the  migration  of  vanadium 
ions  in  the  electrolytes  and  the  oxidation  or  reduction  process  of 
vanadium  ions  at  the  electrodes.  Thereby,  it  is  of  great  significance 
to  understand  the  kinetic  processes  of  redox  reactions  at  the  elec¬ 
trodes  of  VRB  for  developing  electrode  materials  and  optimizing 
the  VRB  structure. 

In  comparison  with  the  negative  V2+/V3+  couple,  the  redox  re¬ 
action  for  the  positive  V0f/V02+  couple  is  more  complicated,  with 
a  multistep  reaction  process  involving  an  oxygen  transfer  reaction 
before  or  after  an  electron-transfer  step  [1].  Since  the  concept  of 
VRB  was  proposed  by  Skyllas-Kazacos  et  al.  in  the  1980s  [4,5],  only 
a  limited  number  of  investigations  on  the  redox  reaction  mecha¬ 
nism  of  the  V0j/V02+  couple  have  been  available  [6—11  ].  In  a  study 
of  the  redox  reaction  mechanism  of  the  V0j/V02+  couple  at  a 
rotating  disk  graphite  electrode  in  acidic  aqueous  solutions,  Gattrell 
et  al.  6]  observed  unusually  high  Tafel  slopes  for  the  reduction  of  V 
(V),  and  suggested  that  the  oxidation  reaction  of  V02+  follows  a 
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multistep  electrochemical-chemical-chemical  (ECC)  mechanism, 
while  the  reduction  reaction  mechanism  of  VOf  changes  from  a 
multistep  chemical-electrochemical-chemical  (CEC)  mechanism  at 
low  cathodic  overpotentials  to  a  multistep  electrochemical- 
chemical-chemical  (ECC)  mechanism  at  high  cathodic  over¬ 
potentials.  However,  the  unusually  high  Tafel  slopes  (350- 
450  mV  decade-1)  observed  for  the  reduction  of  VOf  at  higher 
overpotentials  could  not  be  well  explained  by  the  model.  To  un¬ 
derstand  the  unusual  Tafel  slope,  Gattrel  et  al.  [7]  further  investi¬ 
gated  the  electrochemical  reduction  of  VOj  at  high  overpotentials 
by  polarization  and  double  layer  measurements,  and  established  a 
model  involving  electron  transfer  through  an  adsorbed  layer  of 
vanadium  species  to  explain  the  change  of  the  reduction  reaction 
mechanism  of  VOj  [7].  The  experimental  data  for  the  redox  re¬ 
actions  of  the  V0f/V02+  couple  on  a  rotating  graphite  disk  elec¬ 
trode  have  fitted  well  with  the  computational  simulation  results 
based  on  the  above  model  to  some  extent  [6,7].  Additionally,  Oriji 
et  al.  [11  ]  reported  that  the  redox  processes  of  the  V0f/V02+  couple 
at  a  carbon  plate  electrode  in  both  directions  are  not  simple  and  a 
probable  reversible  chemical  step  with  the  breaking/formation  of 
the  V— O  chemical  bond  follows  an  electron  transfer  step  of  the 
V0j/V02+  reaction  in  a  V  (IV)  or  V  (V)  solution.  Thereby,  more  in¬ 
vestigations  are  still  required  to  understand  the  mechanism  and 
kinetic  processes  of  the  redox  reactions  of  the  V0f/V02+  couple  on 
carbon  electrodes. 

In  this  paper,  the  oxidation  reaction  of  V02+  has  been  investi¬ 
gated  by  steady-state  potentiodynamic  polarization  on  rotating 
graphite  and  glassy  carbon  disk  electrodes  in  sulfuric  acid  solutions 
with  various  pH  and  vanadium  concentrations  in  an  attempt  to 
establish  its  reaction  mechanism. 

2.  Experimental 

2.1.  Preparation  and  chemical  analysis  of  electrolyte  solutions 

All  chemicals  used  in  the  present  experiments  were  analyti¬ 
cally  pure  agents  and  all  solutions  were  prepared  with  de-ionized 
water.  V  (IV)  acidic  solutions  were  produced  by  dissolving 
V0S04.nH20  (n  =  2.82,  according  to  the  chemical  precipitation 
method)  in  sulfuric  acid  solutions  and  then  diluted  to  produce 
the  solutions  with  various  pH  and  vanadium  concentrations.  In 
addition,  the  measurements  of  pH  and  viscosity  of  the  electro¬ 
lytes  were  performed  with  a  pHS-25  pH  meter  (Shanghai  Preci¬ 
sion  Instrument  Co.  Ltd,  China)  and  a  NDJ-5S  digital  torsion 
viscosimeter  (Shanghai  Ni  Run  Intelligent  Technology  Co.  Ltd, 
China),  respectively. 


Instrument  (Shanghai  Chenhua  Instrument  Co.  Ltd,  China). 
Additionally,  some  potentiodynamic  polarization  curves  were 
also  obtained  on  a  spectroscopically  pure  graphite  electrode 
(0.283  cm2,  Sinosteel  Shanghai  Advanced  Graphite  Material  Co. 
Ltd,  China)  at  the  scan  rate  of  0.5  mV  s-1  in  static  V  (IV)  acidic 
solutions  at  various  temperatures  controlled  with  a  water  bath. 
The  electrode  was  ground  with  down  to  2000  grit  silicon  carbide 
papers,  followed  by  polishing  with  a  flanelle  and  thoroughly 
rinsing  with  alcohol  before  use. 

3.  Results  and  discussion 


3.1.  Effect  of  rotation  rate 

Fig.  1  shows  the  polarization  curves  of  the  rotating  pyrolytic 
graphite  disk  electrode  at  various  rotation  rates  in  0.5  M 
VOS04  +  3.0  M  H2SO4  solution.  The  limiting  diffusion  current  (If) 
increases  with  the  increase  of  rotation  rate,  due  to  the  fact  that  a 
higher  rotation  rate  contributes  a  thinner  diffusion  layer  [12]. 
However,  the  anodic  oxidation  current  changes  little  when  the 
rotation  rate  is  higher  than  1600  rpm,  suggesting  that  a  rotation 
rate  larger  than  1600  rpm  is  favorable  to  study  the  oxidation  re¬ 
action  mechanism  of  V02+.  Thus,  a  rotation  rate  of  2000  rpm  has 
been  selected  for  the  subsequent  measurements  on  the  rotating 
disk  electrode  (RDE). 


3.2.  Reaction  order  measurements 


3.2.1.  Reaction  order  of  V02+ 

Fig.  2a  gives  the  polarization  curves  of  the  pyrolytic  graphite 
electrode  at  the  rotating  rate  of  2000  rpm  in  VOSO4  +  3.0  M  H2S04 
solutions  containing  various  concentrations  of  V02+.  Fig.  2b  shows 
clearly  that  the  shift  of  the  polarization  potential  (E)  with  the  log¬ 
arithmic  current  (lgJ)  follows  a  linear  law  at  potentials  from  0.90  to 
1.02  V,  suggesting  that  the  oxidation  reaction  of  V02+  is  controlled 
by  the  electrochemical  polarization.  The  reaction  order  of  V02+ 
(mvo2+)  can  be  calculated  at  a  certain  potential  by  the  following 
equation  [13]: 


111  = 


(  aig/  \ 

\dlgCV02+  )  ^ 


(1) 


where  I  is  the  anodic  current  at  a  certain  potential  (E),  CV02+  is  the 
bulk  concentration  of  V02+. 

From  Fig.  2a  and  b,  the  plots  of  lg I  vs.  lgCV02+  at  different  po¬ 
tentials  are  obtained,  as  shown  in  Fig.  2c.  The  change  of  lg/  with 


2.2.  Electrochemical  measurements 

Electrochemical  measurements  were  undertook  with  a 
rotating  disk  system  (Pine  Instruments)  including  basal  plane 
pyrolytic  graphite  (GB)  and  glassy  carbon  (GC)  disc  electrodes 
(0.196  cm2),  an  AFMSRCE  electrode  rotator,  an  AFE6MB  shaft  and 
a  cell,  with  a  platinum  plate  as  the  counter  electrode  and  a 
saturated  calomel  electrode  (SCE)  as  the  reference  electrode.  A 
salt  bridge  was  used  to  eliminate  the  liquid  junction  potential 
between  the  Luggin  capillary  and  working  electrode.  The  elec¬ 
trodes  attached  to  the  rotating  disk  system  were  initially  polished 
with  aqueous  slurries  with  alumina  powders  (average  particle 
diameter:  0.6  pm)  on  a  piece  of  chamois  cloth  and  then  ultra- 
sonically  cleaned  for  5  min  in  alcohol  and  de-ionized  water  in 
turn.  Steady-state  polarization  curves  at  the  scan  rate  of  1  mV  s-1 
were  obtained  with  the  rotating  disk  system  at  room  temperature 
in  V  (IV)  acidic  solutions.  All  the  electrochemical  measurements 
were  performed  by  a  CHI730C  Electrochemical  Analysis 


Fig.  1.  Polarization  curves  at  various  rotating  rates  on  the  pyrolytic  graphite  electrode 
with  the  sweep  rate  of  1  mV  s-1  in  0.5  M  V0S04  +  3.0  M  H2S04  solution. 
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Fig.  2.  Polarization  curves  (a,  b)  at  varying  concentrations  of  V02+  and  the  corresponding  logarithmic  current  as  a  function  of  the  logarithmic  concentration  of  V02+  at  various 
electrode  potentials  (c)  on  the  pyrolytic  graphite  electrode  at  the  rotating  rate  of  2000  rpm  with  the  sweep  rate  of  1  mV  s-1  in  V0S04  +  3.0  M  H2S04  solution. 


lgCV02+  is  approximately  linear,  and  thus  the  values  of  mvo 2+  (the 
slope  of  the  plots)  at  various  potentials  are  estimated,  as  seen  in 
Table  1.  The  mean  value  of  mvo 2+  is  1.92  for  the  oxidation  reaction  of 
V02+  on  the  pyrolytic  graphite  electrode. 

Fig.  3a  and  b  present  the  polarization  curves  of  the  rotating 
glassy  carbon  disk  electrode  with  a  rotating  rate  of  2000  rpm  in 
VOSO4+  3.0  M  H2SO4  solutions  containing  various  concentrations 
of  V02+  and  the  corresponding  plots  of  lg /  vs.  lgCVC)2+  at  different 
potentials,  respectively.  Linear  Tafel  regions  appear  distinctly  on 
the  polarization  curves.  However,  some  non-linear  features  are 
also  observed  on  the  polarization  curves  at  high  overpotentials 
more  than  120  mV,  probably  due  to  the  glassy  carbon  electrode 
with  a  low  activation.  Therefore,  linear  Tafel  regions  at  as  high 
overpotentials  as  possible  are  selected  to  calculate  the  reaction 
orders.  The  mean  reaction  order  of  V02+  (mVC)2+)  is  similarly 
calculated  to  be  1.55,  as  listed  in  Table  2.  Based  on  the  above 
results,  it  is  approximately  thought  that  the  oxidation  reaction 
order  of  V02+  on  carbon  electrodes  is  2.  It  should  also  be  noted 
that  as  compared  to  the  pyrolytic  graphite  electrode  (Fig.  2a),  the 
open  circuit  potential  for  the  glassy  carbon  electrode  shifts 
negatively  by  around  300  mV,  except  for  the  solution  containing 
0.5  M  V02+.  This  may  be  ascribed  to  their  differences  in  the 
structure  and  surface  activation. 


Table  1 

The  reaction  order  of  V02+  during  the  oxidation  reaction  of  V02+  on  the  pyrolytic 
graphite  electrode. 


E/V 

0.92 

0.95 

0.97 

1.00 

1.02 

Mean  value  of  mv02+ 

mvo2+ 

1.88 

1.93 

1.94 

1.94 

1.93 

1.92 

3.2.2.  Reaction  order  of  H+ 

Fig.  4a  presents  the  polarization  curves  for  the  rotating  pyrolytic 
graphite  electrode  in  0.5  M  VOSO4  +  H2S04  solutions  with  various 
values  of  pH.  At  the  potentials  from  0.75  to  0.9  V,  E  is  linearly 
proportional  to  lg/  (Fig.  4b),  suggesting  that  the  electrochemical 
polarization  is  a  rate  determining  process  for  the  oxidation  reaction 
of  V02+.  The  reaction  order  of  H+  (mH+)  can  be  estimated  as  the 
following  equation  13] : 


(2) 


Fig.  4c  gives  the  shift  of  lg/  with  pH  at  various  potentials.  The 
values  of  mH+  can  be  obtained  from  the  slope  of  the  plots  in  Fig.  4c, 
as  listed  in  Table  3.  The  average  value  of  mH+  is  -1.03  for  the 
oxidation  reaction  of  V02+  on  the  pyrolytic  graphite  electrode. 
Similarly,  the  reaction  order  of  H+  on  the  rotating  glassy  carbon 
electrode  has  also  been  investigated,  as  shown  in  Fig.  5  and  Table  4, 
respectively.  Additionally,  some  non-linear  features  on  the  polari¬ 
zation  curves  are  similarly  observed  in  Fig.  5a.  The  average  value  of 
mH+  is  -0.89.  As  a  result,  the  value  of  mH+  can  be  considered  as  -1 
for  the  oxidation  reaction  of  V02+  on  carbon  electrodes. 


3.3.  Theoretical  calculations  of  reaction  order 

The  above  results  indicate  that  the  reaction  order  of  V02+  and 
H+  for  the  oxidation  of  V02+  on  carbon  electrodes  is  2  and  -1, 
respectively.  The  confirmation  of  these  reaction  orders  is  helpful  for 
understanding  the  oxidation  reaction  process  of  V02+.  Before 
establishing  the  reaction  mechanism  of  V02+,  it  is  necessary  to 
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discuss  theoretically  the  reaction  models  available.  Gattrell  et  al.  [6] 
have  proposed  that  the  oxidation  reaction  of  V02+  is  likely  to  follow 
a  three-step  reaction  mechanism  such  as  chemical  (C)-chemical 
(C)-electrochemical  (E),  ECC  and  CEC  mechanism.  In  addition,  Oriji 
et  al.  [11  ]  have  also  suggested  that  the  electron  transfer  step  (E  step) 
is  followed  by  a  certain  chemical  step  (C  step)  for  the  oxidation 
reaction  of  V02+.  Thereby,  in  the  following  discussions  ECC  and  CEC 
mechanisms  will  be  considered. 

3.3.1.  Calculations  based  on  an  ECC  mechanism 

Supposing  the  oxidation  reaction  of  V02+  is  consistent  with  an 
ECC  mechanism,  then  the  electrochemical  and  chemical  steps  are 
expressed  as  follows  6]: 

V02+  •  5H20  ^  V03+  •  5H20  +  e  (3) 

k- 3 

V03+-5H20^V0j-4H20  +  2H+  (4) 

k—4 

In  the  present  calculations,  two  chemical  steps  are  combined  to 
be  one  chemical  step  of  step  (4).  Any  reaction  rate  constants  (/<)  are 
considered  to  be  unconcerned  with  the  electrode  potential  (E)  in 
this  paper,  except  for  special  annotations. 

If  step  (3)  is  rate  determining,  the  reaction  current  density  i  is 
given  by: 


Table  2 

The  reaction  order  of  V02+  during  the  oxidation  reaction  of  V02+  on  the  glassy 
carbon  electrode. 


E/V 

0.87 

0.9 

0.92 

0.95 

0.97 

Mean  value  of  mvo 2+ 

mvo2+ 

1.57 

1.55 

1.55 

1.55 

1.55 

1.55 

i  =  Fl<4  [V03+]  (8) 

when  step  (3)  is  at  equilibrium,  the  same  reaction  rate  in  both 
directions  is  described  by: 


-q-i) 

RT 


FE 


(9) 


where  k-3  is  the  reverse  reaction  rate  constant  of  step  (3). 

According  to  eqs.  (8)  and  (9),  the  reaction  current  density  i  can 
be  expressed  by  eqn.  (10): 


i  =  F 


k4k3 


[vo2+] 


exp 


k-s 

Thus,  the  values  of  mvo: 
9lg  i 


and  mH+  can  be  calculated  as: 


mvo2+  — 


aig[vo: 


,2+1 


=  1 


(10) 


(11) 


i  =  f/<3  |v02+j  exp 


(5) 


mH+ 


al  gi 

aig[H+] 


=  o 


(12) 


where  k3  is  the  forward  reaction  rate  constant  of  step  (3),  [V02+] 
represents  the  bulk  concentration  of  V02+,  (i  is  the  anodic  charge 
transfer  coefficient,  E  is  the  electrode  potential,  the  other  symbols 
such  as  Ff  R  and  T  have  their  usual  meanings. 

Then  the  reaction  orders  of  V02+  (mvo 2+)  and  H+  (mH+)  can  be 
obtained  from  eq.  (5)  as  follows: 


mvo2+  — 


aigi 


0lg[VO 


,2+i 


=  1 


(6) 


mH+ 


al  gi 

aig[H+] 


=  o 


(7) 


On  the  other  hand,  assuming  that  step  (4)  is  rate  determining 
and  step  (3)  is  close  to  equilibrium,  the  reaction  current  density  i  is 
defined  as: 


The  above  calculations  indicate  that  in  case  the  oxidation  reac¬ 
tion  of  V02+  follows  an  ECC  mechanism,  the  reaction  order  of  V02+ 
and  H+  is  1  and  0,  respectively.  It  is  clear  that  the  theoretical  re¬ 
action  orders  are  inconsistent  with  the  experimental  results,  sug¬ 
gesting  that  the  oxidation  of  V02+  may  follow  a  different 
mechanism. 

3.3.2.  Calculations  based  on  a  CEC  mechanism 

According  to  a  CEC  mechanism,  the  oxidation  reaction  of  V02+  is 
as  follows  [6]: 

V02+-5H20  ^  HVOj -4H20  +  2H+  (13) 

k- 13 

HVOJ  -4H20  U  HV02+  -4H20  +  e  (14) 

k~\4 


Fig.  3.  Polarization  curves  at  various  concentrations  of  V02+  (a)  and  the  logarithmic  current  as  a  function  of  the  logarithmic  concentration  of  V02+  at  various  electrode  potentials  (b) 
on  the  glassy  carbon  electrode  at  2000  rpm  with  the  sweep  rate  of  1  mV  s-1  in  VOS04  +  3.0  M  H2S04  solution. 
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HV02+-4H20  ^  V0j-4H20  +  H+  (15) 

k- is 

Firstly,  assuming  that  step  (13)  is  rate  determining,  the  reaction 
current  density  is  illustrated  as: 

i  =  Ffc13  [V02+]  (16) 

Then,  the  values  of  mvo 2+  and  mH+  are  calculated  to  be  1  and  0, 
respectively. 

Secondly,  if  step  (14)  is  rate  determining  and  the  other  steps  are 
close  to  equilibrium,  the  reaction  current  density  is  given  by: 


i  =  Fk 


14 


[HVOjjexp 


'0  FF 

1?Tfe 


(17) 


Because  step  (13)  is  at  equilibrium,  the  forward  reaction  rate  is 
equal  to  the  reverse  one: 


Table  3 

The  reaction  order  of  H+  on  the  pyrolytic  graphite  electrode  during  the  oxidation 
reaction  of  V02+. 


E/V 

0.78  V 

0.8  V 

0.82  V  0.85  V 

0.87  V 

Mean  value  of  mH+ 

mH+ 

-1.03 

-1.02 

-1.02  -1.03 

-1.04 

-1.03 

i  =  Fk15  [HVO^+]  (20) 

The  following  expressions  are  obtained  for  the  reactions  (13) 
and  (14)  at  equilibrium: 


F/<14  [HVOj jexp 


JLfe 

RT 


=  Fk 


-14 


[hvo|+] 


exp 


-(1  -  0) 


RT 


FE 


(21) 


/<13[V02+]  =  /<_13[h+][hVO+]  (18) 

According  to  eqs.  (17)  and  (18),  the  reaction  current  density  is 
expressed  by: 


i  =  F/<14 


k  13  [V02+] 
k- is  [H+] 


exp 


(19) 


From  eq.  (19)  the  values  of  mvo 2+  and  mH+  are  calculated  to  be  1 
and  -1,  respectively. 

Finally,  if  the  reaction  (15)  is  rate  determining,  then  the 
reactions  (13)  and  (14)  can  be  approximately  considered  to  be  at 
equilibrium.  The  reaction  current  density  is  shown  as: 


Thus,  according  to  eqs.  (18),  (20)  and  (21),  the  reaction  current 
density  is  presented  as  below: 


■  15^14^13 

k-!4k-i3 


[vo2+] 

IhT 


exp 


FF' 

RT 


(22) 


It  can  be  estimated  from  eq.  (22)  that  the  values  of  mvo 2+  and 
mH+  are  1  and  -1,  respectively. 

The  above  calculations  show  that  if  the  oxidation  of  V02+  fol¬ 
lows  the  CEC  mechanism  the  theoretic  values  of  myQ 2+  and  mH+  are 
different  from  the  experimental  results,  suggesting  that  the 
oxidation  reaction  of  V02+  does  not  obey  the  CEC  mechanism  as 
well. 


Fig.  4.  Polarization  curves  (a,  b)  at  varying  pH  and  lg I  as  a  function  of  pH  at  various  electrode  potentials  (c)  on  the  pyrolytic  graphite  electrode  at  2000  rpm  with  the  sweep  rate  of 
1  mV  s-1  in  0.5  M  VOS04  +  H2S04  solution. 
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Fig.  5.  Polarization  curves  at  varying  pH  (a)  and  the  plot  of  lg /  vs.  pH  at  various  electrode  potentials  (b)  on  the  glassy  carbon  electrode  at  2000  rpm  with  the  sweep  rate  of  1  mV  s  1 
in  0.5  M  V0S04  +  H2S04  solution. 


3.4.  A  new  consideration  for  the  oxidation  mechanism  of  V02+ 

The  above  experimental  results  suggest  that  the  reaction  order 
of  V02+  and  H+  for  the  oxidation  of  V02+  on  carbon  electrodes  is  2 
and  -1,  respectively.  To  understand  the  reaction  mechanism  it  is 
necessary  to  know  the  electron  transfer  number  (n)  of  the  oxidation 
reaction  of  V02+.  From  Fig.  1,  the  shift  of  Id  at  the  potential  of  1.4  V 
with  the  square  root  of  the  rotation  rate  (w1^2)  was  obtained,  as 
shown  in  Fig.  6.  Fig.  6  shows  clearly  that  the  change  Id  vs.  w1/2  is 
approximately  linear.  The  value  of  n  can  be  determined  from  the 
slope  of  the  plot  by  Levich’s  equation  [12]: 

Id  =  0.62nFAD2/o32+!/-V6wi/2Cvo2+  (23) 

where  A  is  the  area  of  the  electrode,  DV02+  the  diffusion  coefficient 
of  V02+,  v  the  kinematic  viscosity  of  the  electrolyte  and  Cvo 2+  the 
bulk  concentration  of  V02+. 

v  can  be  defined  by  the  following  equation  [14]: 

7]  =  v  x  p  (24) 

where  rj  is  the  kinetic  viscosity  of  the  electrolyte  and  p  is  the  density 
of  the  solution. 

For  the  present  system,  the  values  of  Dvo 2+  and  p  can  be  assumed 
as  1CT6  cm-2  s-1  11]  and  1.2  g  cm-3  [14],  respectively.  The  kinetic 
viscosity  of  the  electrolyte  (■ r /)  is  measured  as  2.85  mPa  s.  Thus,  the 
value  of  n  is  estimated  as  0.81  from  eqs.  (23)  and  (24).  It  can  be 
deduced  that  the  electron  transfer  number  (n)  of  the  oxidation 
reaction  of  V02+  is  approximately  1.  It  has  been  accepted  that  in 
acidic  solutions  with  a  FI2SO4  concentration  less  than  4.0  M,  V02+ 
and  VOj  exist  as  the  octahedrally  coordinated  cation  V02+-5H20 
[14]  and  the  mononuclear  aquadioxo  cation  V0J.4H20  [15,16], 
respectively.  It  is  also  reported  that  dinuclear  oxovanadium  (IV,  V) 
species  containing  the  [V20s]3+core  can  be  easily  formed  from  the 
mixtures  of  the  constituents  V02+  (aq)  and  VOf  (aq)  ions  because 
of  their  favorable  thermodynamic  stability  [17-20].  Flowever,  the 
dinuclear  oxovanadium  (IV,  V)  complex  formed  shows  the  kinetic 
lability  due  to  the  weak  interaction  between  the  V02+  (aq)  and  VOj 


Table  4 

The  reaction  order  of  H+  on  the  glassy  carbon  electrode  during  the  oxidation  reac¬ 
tion  of  V02+. 


E/V 

0.68  V 

0.69  V 

0.7  V 

0.71  V 

0.72  V 

Mean  value  of  mH+ 

mH+ 

-0.91 

-0.90 

-0.89 

-0.87 

-0.85 

-0.89 

(aq)  ions,  and  thus  the  species  V02+  (aq)  and  VOf  (aq)  are  only 
slightly  consumed  during  the  reaction  [17  .  Thus,  the  formation  of 
the  oxovanadium  (IV,  V)  complex  is  out  of  the  consideration  here 
for  simplicity.  Thereby,  based  on  the  experimental  values  of  mvo 2+, 
mH+  and  n,  a  new  CEC  mechanism  involving  two  intermediates 
HV203+-9H20  and  HV02+-4H20  is  proposed  as  follows: 

2[V02+-5H20|  U  HV20^+-9H20  +  H+  (25) 


HV20|+-9H20  ^  HV02+-4H20  +  V02+-5H20  +  e  (26) 

k-24 


HVC>2+  -4H20  =1  VOJ  -4H20  +  H+  (27) 

k-  25 

Because  the  dinuclear  mixed-valence  oxidovanadium  (IV,  V) 
complex  is  considered  to  have  almost  a  linear  V-O-V  bridge  with 
the  terminal  oxo  groups  in  mutually  transpositions  shown  as 
O 

II 

V— O— V  [17-19],  the  intermediates  HV203+-9Fl20  and 

II 

0 

0 

0  II 

HVOj+  4H20  probably  have  the  structural  forms  of  [  v-O-V  ]3+.9H20 

I 

OH 


Fig.  6.  The  shift  of  /d  at  the  potential  of  1.4  V  with  w1/2. 
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lg(i/A.cm‘2) 

Fig.  7.  Polarization  curves  (a)  and  the  shift  of  E  vs.  lgi  (b)  at  various  temperatures  on 
V0S04  +  0.5  M  H2S04  static  solution. 


Ig(i/A.cm'2) 

spectroscopically  pure  graphite  electrode  with  the  sweep  rate  of  0.5  mV  s-1  in  1.0  M 


and 


[O-V]2-.  4H20 
OH 


,  respectively,  which,  however,  have  not  been 


investigated  experimentally  in  the  present  study. 

The  theoretical  values  of  mvo 2+  and  mH+  for  the  oxidation  of 
V02+  based  on  the  reactions  (25)-(27)  are  similarly  calculated  as 
the  above. 

Firstly,  if  step  (25)  is  rate  determining,  the  reaction  current 
density  is  given  by: 


i  =  flc23[v  02+]2  (28) 

From  eq.  (28)  the  values  of  mvo 2+  and  mH+  are  calculated  to  be  2 
and  0,  respectively. 

Next,  in  case  step  (26)  is  rate  determining  and  the  other  steps 
are  approximately  at  equilibrium,  the  reaction  current  density  is 
expressed  by: 


i  =  Fk24[HV20|+]exp 


(29) 


For  step  (25)  at  equilibrium,  the  following  expression  can  be 
obtained: 


/<23  [V02+] 2  =  k_ 23  [h+]  [HV202+]  (30) 

Introducing  eq.  (30)  into  eq.  (29),  eq.  (29)  is  changed  to: 


i  =  Fk24 


k23  [V02+] 

k- 23  [H+] 


2 


exp 


(31) 


From  eq.  (31 )  the  values  of  mvo2+  and  mH+  can  be  estimated  to  be 
2  and  -1,  respectively. 

Finally,  if  step  (27)  is  rate  determining,  the  steps  (25)  and  (26) 
can  be  thought  to  be  at  equilibrium.  Then,  the  reaction  current 
density  is  illustrated  by: 

i  =  Fk25[  HV02+]  (32) 


Table  5 

The  charge  transfer  coefficients  (0)  on  the  spectroscopically  pure  graphite  electrode 
at  various  temperatures. 


T/°C 

10 

20 

30 

40 

Mean  value 

Tafel  slope  (0.68-0.78  V)/(V  dec”1’ 
0 

)  0.111 
0.531 

0.116 

0.509 

0.107 

0.551 

0.114 

0.518 

0.53 

The  same  reaction  rate  in  the  forward  and  reverse  directions  for 
both  the  steps  (25)  and  (26)  can  be  described  as: 


Fk24[HV20!+] 


exp 


—FE 

RT 


f/<_24[hvo2+] 

x  |v02+j  exp 


-(1  -  0) 


RT 


FE 


(33) 


'<23  [V02+] 2  =  k_ 23  [H+]  [HV202+]  (34) 

Thus,  according  to  eqs.  (32)— (34),  the  reaction  current  density 
can  be  given  by: 


I  =  pk 25^24^23 
k~ 24^—23 


[vo2+] 

~WT 


exp 


FE 

RT 


(35) 


It  can  be  estimated  from  eq.  (35)  that  the  values  of  mvo2+  and 
mH+  are  1  and  -1,  respectively. 

The  above  calculations  indicate  that  when  the  step  (26)  (E  step) 
is  rate  determining,  the  theoretical  values  of  mvo 2+  and  mH+  comply 
well  with  the  experimental  results.  Thereby,  it  is  reasonable  to 
consider  that  the  oxidation  reaction  of  V02+  may  follow  the  new 
CEC  mechanism  with  an  electrochemical  reaction  as  the  rate- 
limited  process.  In  contrast  with  the  mechanisms  proposed  by 
Gattrel  et  al.  [6,7],  the  new  CEC  mechanism  involves  a  possible 
combination  between  V02+  ions  during  the  oxidation  reaction  of 
V02+.  Moreover,  the  reaction  orders  of  V02+  and  H+  for  the 
oxidation  reaction  of  V02+  are  2  and  -1,  respectively.  In  this  case, 
the  reaction  current  density  of  V02+  shown  in  eq.  (31)  can  be  also 
presented  by 


Table  6 

The  formal  potential  (£°  )  and  the  rate  constants  of  the  vanadium  redox  reaction  (ks) 
at  different  temperatures  calculated  from  the  cyclic  voltammograms  for  the  VO2  / 
V02+  redox  reaction  on  the  spectroscopically  pure  graphite  electrode. 


T/°C 

E°'IV 

ksl(  10  3  cm  s  3) 

10 

0.940 

1.91 

20 

0.940 

2.13 

30 

0.935 

3.26 

40 

0.945 

5.03 

50 

0.935 

5.11 
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i  =  Fk 


[V02+]2 

~wr 


exp 


(36) 


where  k  is  the  reaction  rate  constant.  When  the  electrochemical 
reaction  rate  constant  at  standard  state  (/<s)  is  taken  as  a  substitute 
for  /<,  eq.  (36)  can  be  described  as  [12]: 


i  =  Fks 


[vo2+]2 

~wr 


exp 


(37) 


where  E°  is  the  formal  potential  of  the  V02+/V0f  redox  reaction. 


3.5.  An  oxidation  reaction  kinetic  equation  of  V02+  on  the 
spectroscopically  pure  graphite  electrode 

The  above  results  indicate  that  the  proposed  CEC  model  with  the 
electrochemical  reaction  as  the  rate-limiting  step  is  reasonable  to 
be  used  to  describe  the  oxidation  reaction  of  V02+.  To  establish  the 
real  reaction  kinetics  of  V02+,  some  parameters  shown  in  eq.  (37) 
should  be  known.  Fig.  7a  shows  the  polarization  curves  at  various 
temperatures  on  the  spectroscopically  pure  graphite  electrode  in 
1.0  M  V0S04  +  0.5  M  H2S04  static  solutions.  Fig.  7b  presents  the 
corresponding  plots  of  E  vs.  lgi  with  the  slope  as  Tafel  slope  (9E/ 
dlgi).  Furthermore,  the  parameter  in  eq.  (37)  can  be  defined  as  the 
following  equation  [12]: 


dE  _  2.303 RT 
01  gi  -  |3nF 


(38) 


Thus,  the  values  of  Tafel  slope  and  j3  at  different  temperatures 
can  be  obtained,  as  listed  in  Table  5.  The  mean  value  of  /3  is  obtained 
to  be  0.53. 

With  the  cyclic  voltammograms  on  the  spectroscopically 
pure  graphite  electrode  in  1.2  M  VOSO4  +  3.0  M  H2S04  static 
solutions  at  various  temperatures  examined,  the  formal  poten¬ 
tial  (F° )  and  the  rate  constants  of  the  vanadium  redox  reaction 
at  standard  state  (ks)  at  various  temperatures  are  estimated,  as 
listed  in  Table  6,  and  the  corresponding  calculations  have  been 
shown  in  detail  in  another  unpublished  paper  submitted.  The 
mean  value  of  E°  is  0.94  V.  Based  on  Table  6,  the  change  of  In  ks 
with  -1/T  is  obtained,  as  shown  in  Fig.  8.  The  reaction  rate 
constant  ks  at  standard  state  can  be  given  by  the  Arrhenius 
equation  [12]: 


ks  =  A0  exp 


(39) 


where  A G6  is  the  reaction  activation  energy  at  standard  state,  Ao 
the  temperature-independent  factor.  Thus,  from  Fig.  8,  the  values 
of  A Gd  and  A0  can  be  estimated  as  21.774  kj  mol-1  and 
18.202  cm  s-1,  respectively.  Introducing  the  values  of  these  pa¬ 
rameters  into  eq.  (37),  the  kinetic  equation  for  the  oxidation 
reaction  of  V02+  on  the  spectroscopically  pure  graphite  electrode 
can  be  given  as: 


i  =  1.756  x  lO^^C^1  exp 


69820\ 


RT 


0.53  FE\ 


exp 


(40) 


where  E  is  the  electrode  potential  versus  SCE,  i  is  the  anodic  current 
density  (A  cm-2),  Cvo 2+  and  CH+  are  the  bulk  concentrations  of  V02+ 
and  H+  ions  in  the  electrolyte  solutions  (mol  cm-3),  respectively. 

According  to  the  kinetic  eq.  (40),  the  change  of  lgi  with  E  can 
be  obtained.  Fig.  9a  presents  a  comparison  between  the  theo¬ 
retical  and  experimental  polarization  curves  at  10  °C  in  1.0  M 
VOSO4  +  0.5  M  H2SO4  static  solution.  It  can  be  seen  that,  except 
for  the  potentials  higher  than  0.78  V,  the  real  kinetic  equation  fits 
the  experimental  data  quite  well.  Since  the  kinetic  equation  is 
just  usable  at  the  potentials  where  the  electrochemical  polari¬ 
zation  is  rate  determining,  these  are  the  very  expected  results. 
Similar  results  have  also  been  obtained  at  25  °C  in  1.5  M 
VOSO4  +  3.0  M  H2SO4  static  solution,  as  shown  in  Fig.  9b, 
respectively.  Therefore,  the  kinetic  equation  of  eq.  (40)  can  be 


Fig.  9.  The  kinetic  equation  fitted  to  a  set  of  polarization  plots  on  the  spectroscopically  pure  graphite  electrode  in  1.0  M  V0S04  +  0.5  M  H2S04  static  solution  at  10  °C  (a)  and  1.5  M 
VOSO4  +  3.0  M  H2S04  static  solution  at  25  °C  (b). 
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well  used  to  predict  the  polarization  behavior  of  the  spectro¬ 
scopically  pure  graphite  electrode  over  a  range  of  temperatures 
and  concentrations  in  V  (IV)  acidic  solutions. 

4.  Conclusions 

The  oxidation  reaction  mechanism  of  V02+  on  a  rotating  disk 
electrode  such  as  basal  plane  pyrolytic  graphite  (GB)  and  glassy 
carbon  (GC)  disk  electrodes  has  been  examined  by  potentiody- 
namic  polarization  in  sulfuric  acid  solutions  with  various  pH  and 
vanadium  concentrations,  respectively.  Based  on  the  polarization 
results,  the  reaction  orders  of  V02+  (mvo 2+)  and  H+  (mH+)  for  the 
oxidation  reaction  of  V02+  on  carbon  electrodes  have  been 
calculated  to  be  2  and  -1,  respectively,  which  are  inconsistent 
with  the  theoretical  results  based  on  the  possible  reaction 
mechanisms  including  chemical  (C)-chemical  (C)-electro- 
chemical  (E),  ECC  and  CEC  mechanism  available  in  the  literature. 
The  electron  transfer  number  (n)  of  the  oxidation  reaction  of 
V02+  is  calculated  to  be  approximately  1.  Based  on  the  experi¬ 
mental  values  of  mvo 2+,  mH+  and  n,  a  new  CEC  oxidation  reaction 
mechanism  shown  in  eqs.  (23)— (25)  is  thus  proposed  to  repre¬ 
sent  the  oxidation  of  V02+.  The  theoretic  reaction  orders  of  V02+ 
and  H+  based  on  the  new  mechanism  are  in  good  accord  with 
the  experimental  results  when  the  electrochemical  polarization 
is  rate  determining.  Furthermore,  a  corresponding  kinetic 
equation  has  been  built  up  for  the  oxidation  reaction  of  V02+  on 
a  spectroscopically  pure  graphite  electrode  as  shown  in  eq.  (40), 
and  can  be  well  used  to  predict  the  polarization  behavior  over  a 
range  of  temperatures  and  concentrations  in  V(IV)  acidic  solu¬ 
tions.  Though  the  new  reaction  model  can  be  well  applied  to  a 
certain  extent,  more  extensive  studies  are  still  required  to  un¬ 
derstand  the  new  mechanism,  including  the  oxidation  reaction 
on  different  carbon  electrodes,  the  reaction  intermediates 
involved,  etc. 
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